Prom classical to quantum Kagome antiferromagnet in a magnetic field 
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We study the magnetic properties of the Kagome antiferromagnet going from the classical limit 
to the deep quantum regime of spin 1/2 systems. In all the cases there are special values for the 
magnetization, 1/3 in particular, in which a singular behavior is observed to occur in both the 
classical and quantum cases. We show clear evidence for a magnetization plateau for all S, in a wide 
range of XX Z anisotropics and for the occurrence of quantum order by disorder effects. 
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Spin systems in low dimensions have become an in- 
tense area of research in the last few years due to their 
possible relevance to high- Tc cuprates but also in their 
own right due to their very special properties Q. One 
of the issues that attracted attention, both from exper- 
imentalists and theorists, is the appearance of plateaux 
in the magnetization curves in both ID and 2D 

1^-^ quantum spin systems. Another has been the study 
of strongly frustrated systems with highly degenerate 
classical ground state structure (lO| giving disordered, 
spin liquid behavior at low temperature. The spin liq- 
uid behaviour can be removed, or partial removed by 
degeneracy lifting, through thermal, or quantum fluc- 
tuations; i:he phenomenon known as order by disorder 
(OBD) Examples are the so-called Ji — J2 model 

on a square lattice and the Kagome Heisenberg an- 
tiferromagnet (AF), where thermal OBD leads to a ne- 
matic spin liquid state ]T^-p^. The situation is more 
complicated when one considers quantum spins defined 
on such lattices, and the role of quantum fluctuation in 
the selection of some particular ground states (GS's). In 
the particular case of the Ka gom e lattice , state of the 
art numerical computations |Iq] seem to indicate that no 
nematic order is observed and that the GS is a rather 
complicated spin liquid like state, with a spin gap but 
with many low lying spinless states. 

Recent studies of the Kagome AF in the presence of a 
magnetic field have shown that the value of 1/3 for the 
normalized magnetization is a special point bringing into 
play both OBD effects and the presence of a magnetic 
plateau, an effect also present for the triangular lattice 
In the classical case [^,^, spin wave analysis 
and Monte Carlo simulation show the order by disorder 
selection of a collinear spin-liquid state together with a 
magnetic anomaly and a pseudo magnetic plateau. In the 
quantum case, numerical evidence from the exact diago- 
nalization of small S = 1/2 clusters show the existence 
of a magnetic plateau [ p2[ . Such behavior has already 
been observed in the Ji — J2 square lattice, where similar 
phenomena, based on a collinear order of the spins for 



particular values of the magnetization, occurs | |20[ | . 

A natural question one can pose is therefore: is there 
a connection between the existence of soft modes in 
the classical case and the appearance of magnetization 
plateaux in the quantum system ? Further; is the mag- 
netic plateau accompanied by quantum order by disorder 
selection of a reduced manifold of states ? In this pa- 
per we have addressed these questions using a standard 
spin wave analysis, Lanczos diagonalization of finite clus- 
ters and a 2D Jordan- Wigner approach. We study the 
XXZ AF on the Kagome lattice, where in general no 
soft modes are present. Our results show the persistence 
of the plateau at 1/3 beyond the symmetric XXX point, 
and even the existence of a plateau for classical spins, in 
the easy axis case. This suggests that the existence of the 
soft modes is not a necessary condition for the existence 
of a plateau in the quantum case, although it does seem 
to extend the plateau region away from the Ising limit of 
XXZ AF. The quantum ground state does indeed have a 
large projection onto a collinear state indicating a quan- 
tum order by disorder effect. 

We define the XXZ AF Hamiltonian 
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where Si — {Sf, Sf , Sf) is a vector of unit length, J > 
is the exchange constant and < A < 00 the "XXZ" 
anisotropy parameter. Specific values, A = 1 and A = 
define the XXX and XX models respectively. 

Considering first the XXX case and defining Sp = 
2 3 ^ total moment per triangular unit, the 
Hamiltonian for N spins can be written as 



H = -NJS^ 



2N/3 

E 

p=i 



2 p 



(2) 



Minimizing with respect to Sp gives Sp — h/{2J) and a 
net GS magnetization M = h/{6JS), for \h\ < he — 6JS. 
Above this field the spins are ferromagnetically aligned. 



1 



In zero field this condition leads to a rigid coplanar spin 
configuration for each triangular plaquette, but the open 
nature of the Kagome lattice allows freedom in the orien- 
tation of the triangular spin planes and a highly degen- 
erate continuous GS manifold. Spin wave analysis shows 
that a branch of modes for fluctuations out of a discrete, 
but macroscopically large subset of collectively coplanar 
ground states become soft, leading to their thermal se- 
lection |l|-|l6|,||. 

For non-zero magnetic field, a generic classical ground 
state for a single plaquette is not necessarily coplanar, 
in contrast to the zero field case [Q. Selection of a col- 
lective coplanar state therefore requires selection at the 
level of each triangle, as well as between triangles. 

This fact weakens, but does not destroy the coplanar 
selection For small field Zhitomirsky has proposed 
that a discrete set of coplanar states are selected with, 
on each triangle one spin, A, antiparallel to the field and 
spins B and C straddling the field direction in a "F" 
configuration. As in zero field, these states have soft 
modes associated with out of plane spin fluctuations and 
corresponding line defects in real space, allowing rotation 
of spins B and C about A. 

As h ^ hc/S the two canted spins close up onto the 
field direction, forming coUinear "up-up-down" (UUD) 
spin configurations on each triangle, with magnetization 
M = 1/3. This special point has zero modes for in plane, 
as well as out of plane spin fluctuations |^,|l[ . The large 
number of soft modes leads to strong thermal selection 
of the UUD spin hquid @,|l[. 

To describe the fluctuations over the coplanar states 
near /ic/3, we choose a coordinate system oriented with 
respect to local spin directions Q. At each site we 
choose right handed axes in spin space with the "z" axis 
parallel to the spin and all "x" axes mutually parallel, 
that is, perpendicular to the spin plane. The spin orienta- 
tions are therefore parameterized hy Si = (ef , ej', 1 — a^), 
with ai determined from \Si\ = 1 and the quadratic part 
of the Hamiltonian reads 



1/2 y: m,^^^ 



i,j t J J 



(3) 



The only non-zero elements of the matrices and 
are the diagonal and nearest-neighbors terms. In partic- 
ular, entries are Mf^ = 1 and Mf^ = 1/2 if i and j 
are nearest-neighbors. While turns out to be inde- 
pendent of h and the particular GS we are considering, 
this is not the case for M^. 

The analysis of the "x" sector is given in reference 
[0. Even for disordered states the fundamental excita- 
tions are propagating modes, rather than the localized 
line defects that provide an alternative description of the 
soft fluctuations |Q . One can diagonalize @) in Fourier 
space finding hidden branches of magnon-like excitations, 
as far as magnetic scattering techniques are concerned: 
an optical branch, an acoustic branch and a soft branch 
extending over the entire Brillouin zone. This branch 



is responsible for the selection of coplanar states in zero 
field and is present for all field values [|l^,^. At the 
point h/hc — 1/3 and for fluctuations around an UUD 
state the "y" sector also acquires a soft branch, indepen- 
dently of the GS. This fact plays in favor of a selection 
of the UUD states and leads to the UUD "spin-liquid" 
predicted and observed in reference pG]. 




FIG. 1. Zero temperature magnetization curve, for classical 
spins. A — 1.5. 




2 X = h/(JS) 

FIG. 2. Schematic classical phase diagram in the h, A 
plane. 

The singular entropy of fluctuations driving the UUD 
spin liquid also leads to a magnetic anomaly and the be- 
ginnings of a magnetic plateau. Close to h — hc/3 the 
growing entropic contribution to the free energy pulls 
the minimum away from the ground state energy and 
towards the UUD point. This pulls the magnetization 
above the linear curve M — h/hc as one approaches hc/i 
from below and pulls it below as one approaches from 
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above. However the spin wave analysis gives no plateau, 
rather it gives a continuous relaxation away from linear 
behaviour and a point of inflection at h — hc/S. A re- 
maining open question is whether nonlinear corrections 
call give a discontinuity and a finite plateau, as proposed 

in g. 

For A > 1 and in zero field the ground states re- 
main highly degenerate and there are soft modes for out 
of plane fluctuations about collectively coplanar states 
with moment M = (A - 1)/3(A -f- 1) Adding 
a small field along the z axis does not lift the degen- 
eracy or remove the soft modes. The magnetization 
increases linearly until M = 1/3 at h — 2JS, inde- 
pendently of A. For higher fields, there is a plateau 
in the magnetization, at M = 1/3, continuing out to 
h = (2 A - 1 + \/4A2 + 4A - 7)JS where it starts to 
grow again, until reaching saturation at h = (2 + 4A) JS*. 
The plateau occurs because, in this field region the field 
and exchange contributions to the Hamiltonian both vary 
quadratically with deformation away from an UUD con- 
figuration, leading to a linear instability. As a result 
the macroscopic, but discrete set of UUD states is se- 
lected energetically rather than entropically and the field 
must exceed h = 2JS by a finite amount, before the 
ground state manifold leaves these configurations. The 
in-plane soft modes at M = 1/3 are however a property 
of the XXX point and are removed for A > 1. The 
lower boundary of the plateau tends to Af = as one ap- 
proaches the Ising limit, A — > cxj, and there is a zero tem- 
perature plateau extending over the whole of the range 
h < he- The zero temperature magnetization for the clas- 
sical case, is shown in Fig. 0, for A = 1.5. The data were 
obtained from zero temperature Monte Carlo, decreasing 
the field the starting value h = he- The plateau region is 
clearly visible. For A < 1 and zero field, the out of plane 
fluctuations acquire a finite stiffness |Q. In finite field 
the ground states are formed by canting the spins out of 
the plane perpendicular to the field direction. There are 
in general no soft modes and M = 1/3 no longer has any 
special significance. 

Classically, there is therefore an extended region in the 
A, h phase diagram at low temperature where plateau 
phenomena can be observed, extending from the point 
A = 1, h = 2JS in a wedge, out to the line A = oo, h < 
he, as is shown schematically in Fig. A natural ques- 
tion to ask is therefore, do quantum fluctuations extend 
this region, leading to a plateau over a larger region of 
parameter space? 

In trying to understand this problem, we present here 
a brief description of quantum fluctuations over copla- 
nar states for arbitrary values of M. Within a Holstein- 
Primakoff description for the spin variables, where the 
quantization axis is chosen according to the classical ori- 
entation of the spin described above, the XXX Hamil- 
tonian can be written: 



where H2 is a quadratic Hamiltonian of magnon creation 
and annihilation operators. The 0{l/^/S) part contains 
higher order terms in these operators. By expressing the 
creation and annihilation operators in terms of the cor- 
responding canonical variable Q and P, H2 reads pq] : 



H.2 - 1/2 ^ \MlP^P, + M^Q^Q,] 



(5) 



where and are the matrices for classical fluctua- 
tions in the y and x sectors discussed above. There are 
therefore N/?) Q coordinates absent in H2 with the cor- 
responding canonical variables P commuting with H2. 
This is the quantum analog of the flat branch. It can 
be explicitly shown to exist for the ordered coplanar GS 
with (7 = ordering vector, using a three color boson pic- 
ture in momentum space. In this case, by an appropriate 
canonical transformation, H2 can be written as: 



H. = 



^Y.[p^{k)P,i-k)+uj^ik)X^ik)X^i-k)] , (6) 



k,fi 



where /i is the color index and uj3{k) = V k. We see 
then the presence of the flat mode in the quantum case, 
at this order in 1/y/S, as a free particle like spectrum. 

Of particular interest are the UUD states, where N/3 
P coordinates are also absent from the Hamiltonian at 
this order. It is obvious that higher order terms of the 
Hamiltonian are crucial for understanding the spectrum 
of these excitations. In |2^ it was shown how the zero- 
point fluctuations of the fast variables (here i — 1,2) 
introduce (via the higher order terms) an effective po- 
tential for the soft modes, apart from lifting the classical 



, processes like the tunnel- 
iscrete set of local minima 



S 



H = Ho + -{H2 + 0{\/^/S)) 



(4) 



GS degeneracy. Following 
ing effect between the, now c 
can give rise to a finite, and even gaped stiffness of the 
soft modes. What we should point out in this example, 
is that the case \M\ = 1/3 corresponds to a different sit- 
uation, where the effective dynamics of the soft modes 
introduced by the interactions with the normal oscilla- 
tors will be radically different. We argue that, because 
of the absence of this huge number of canonical variables 
in H2 for the UUD states, quantum fluctuations will in- 
deed play in favor of a selection of these states, and to 
their stability as the magnetic field is changed away from 
hc/3. Such a selection, for a magnetic UUD state, is also 
present for the triangular lattice, due to zero-point fluc- 
tuations over an ordered classical ground state Jl9| , p8t . 
Following this idea one could follow a standard semiclas- 
sical expansion, and argue that if quantum fluctuations 
select a "Y" coplanar configuration for M < 1/3 and a 
quasi-coUinear configuration for Af > 1/3 as do thermal 
fluctuations in the classical case , the value of the field 
at which M — 1/3 is reached from above and below will 
be renormalized at the quantum level. This scenario gives 
rise to a plateau at M = 1/3 even for A = 1. Indeed, tak- 
ing for example a q — structure, a straightforward nu- 
merical integration of the free bosonic branches involved 
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either in the ""K" or quasi-coUinear Kagome vacua, yields 
a finite plateau width. While this estimation is based on 
a particular ground-state {q — in this case), it is rea- 
sonable to assume a similar behavior for any coplanar 
ground-state. A more subtle point concerns the fact that 
this approach has to be taken with some care for the 
Kagome model, because the oscillator frequencies have 
a non-analytic behavior at h = 2JS. This can give rise 
to divergencies in the coefficients of the formal power se- 
ries and produce a non-analytic dependence on S of the 
width of the plateau. 

Motivated by these semiclassical arguments in favor 
of a plateau at 1/3, we further analyze this issue using 
Lanczos diagonalization and a Jordan- Wigner approach, 
which are more appropriate in the deep quantum case 
pof . Let us first focus attention on the exact GS's ob- 
tained from the diagonalization of the clusters with peri- 
odic boundaries schematized in Fig. 3. Using a Lanczos 
algorithm applied through all subspaces, we ob- 
serve a plateau at M — 1/3 as a function of magnetic 
field. Results for = 18 and all 24-spin clusters are 
shown in Fig. 4 for both XXX and XXZ with A > 1. 








FIG. 3. Schematic view of clusters considered in exact di- 
agonalizations. (a) 18 spins; (b), (c) and (d) 24 spins; (e) 
27 spins. Dotted lines indicate PBC, holding throughout all 
boundaries. 




FIG. 4. Magnetization curves for XXZ clusters (A = 2.5) . 
Dashed, solid, short dashed, and dotted lines denote respec- 
tively the results obtained for clusters (a), (b), (c) and (d). 
The inset displays the corresponding results for isotropic ex- 
changes (A = 1) . 

To check whether an UUD order actually occurs, we 
measured three-point coplanar correlation functions of 
the form Ce{i, j, k) = J2<i.j,k> ( "-j" ~ ^k) ) the 
M =1/3 GS. Here, < i, j, k > denotes a cyclic permuta- 
tion around an elementary triangle whereas, — a^a^ 
is the usual spin occupation number and fi^ stands for the 
local density fields along the 6'-axes forming the coplanar 
"y" configuration. Specifically, after some elementary 
steps it can be readily shown that 



<i,j,k> 



{ [ ~ sin^ Q erf (1 cos Q a[) 



X (l + cos0a|)] (1-nD) 



(7) 



It turns out that Ce increases monotonically with A and 
near the maximum depends very little either on the clus- 
ter size, its particular shape, or the plaquette location. A 
large but unsaturated value of this parameter at 6 — , 
as is shown in Fig. 5, is consistent with the idea of a 
collinear UUD spin liquid allowing a finite probability 
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for configurations with defects. 





FIG. 5. Triangular correlation functions measuring varia- 
tions around the coplanar "Y" configurations at M = f/3 in 
cluster (d). Solid, dashed and dotted lines refers to A = 2.5, 1, 
and respectively. For comparison, the inset shows the cor- 
responding situation a, A'l — 0, favoring a classical behavior. 

One can follow this idea further by testing if the sys- 
tem selects a single state or microscopic number of states 
making up the UUD spin liquid. By sorting the proba- 
bilities of the ground state components in real space, we 
find that the latter exhibit a sharply peaked distribution, 
at the top of which there are a number of UUD states. 
This is shown schematically in Fig. 6 for cluster (e). The 
weight distribution becomes even sharper for A > 1 and 
broadens steadily upon decreasing the anisotropy. Sim- 
ilar results were observed in the remaining clusters dis- 
played in Fig. 3 . This ensemble of results therefore 
gives a strong indication that there is quantum order by 
disorder selection of the UUD spin liquid over the plateau 
region for A = 1. The narrowing for A > 1 is evidence 
for the removal of non-coUinear fluctuations away from 
the symmetric XXX point, as one might expect from 
the energetic arguments presented earlier. 



FIG. 6. Ground state of cluster (e) for M = 1/3 and 
A — 2.5 . Real space states are sorted by increasing weights 
and normalized to their maximum value (up-up-down config- 
urations). The inset exhibits the ground state components 
for A = 1 . Only weights above 10"'^ are shown. 

As a final step, we now study our system by means of a 
generalization of the well known Jordan- Wigner transfor- 
mation in ID, which was originally proposed in and 
further developed in ||3^-|35|] . It consists of a mapping of 
the original spin variables into spinless fermion variables 
Q via S~ = CiexpiJ^j^^CjCj (irg{Tj - Ti), where n = 

Xi + iyi is the complex coordinate and Sf — c\ci — 1/2. 
The advantage of this transformation is that it imposes 
exactly the single occupancy constraint, but at the price 
of introducing a minimal coupling with Chern-Simons 
gauge fields. One usually proceeds by a mean field de- 
coupling to end up with a problem of fermions on the 
given lattice in a properly chosen magnetic background, 
which can be solved using standard techniques. In this 
language, an external magnetic field applied to the spin 
system enters as a chemical potential for the spinless 
fermions. In the present case and mainly motivated by 
the tendency of the AF Kagome lattice to disorder , we 
have chosen a homogeneous mean field, that is a uniform 
fermion density cjcj = ^(1 +M), Vj . We restricted 
the analysis to the case A = where this calculation 
scheme becomes most reliable. The resulting magnetiza- 
tion curve is shown in Fig. 7. A clear plateau again shows 
up at M ~ 1/3. Two other plateaux also appeared at 
M — and M = 2/3 which could be an artifact of the 
mean field approximation and more work is needed to be 
conclusive at this point. 
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FIG. 7. Mean field magnetization curves for XX lattices, 
using homogeneous fermion densities. For comparison, the 
inset shows the exact magnetization curve of cluster (c) using 
A = 0. 

It is interesting to recall that the 1/3 plateau is ex- 
pected to disappear for low enough XX Z anisotropy in 
the triangular lattice ||^ , while our mean-field results sug- 
gest that it persists even in the XX limit of the Kagome 
lattice. This result is again consistent with the idea that 
quantum fluctuations stabilize a plateau over an exten- 
sive range of parameter space. The Lanczos magnetiza- 
tion curve at A = yields no conclusive evidence on 
this issue, though there is evidence that the sharp weight 
structure of the GS components disappears. 

To summarize, we have studied the magnetization be- 
havior of both classical and quantum Kagome AF using a 
variety of methods. In all cases our results suggest that at 
least for A > 1 , there is selection of an UUD state which 
is ultimately responsible for the emergence of a magne- 
tization plateau at M — 1/3. An additional feature, 
particular to the Kagome lattice is that the UUD state 
is not a Neel ordered state, rather it is a spin liquid. Our 
evidence for this is that we find a number of states with 
comparable weights. It seems therefore that the strong 
fluctuations associated with the the spin liquid state do 



not destroy the plateau; rather, the two effects co-habit. 
Further investigation of these issues is underway. At the 
classical level a plateau occurs in a region of (A, h) phase 
space, terminating at the point (A = 1, h = 2JS), where 
a pseudo plateau is driven by the entropic contribution 
to the free energy. The results of this paper allow us to 
speculate that quantum fluctuations stabilize a plateau 
over a more extensive region of parameter space. The 
robustness of this result appears to be associated with 
a clean separation of the GS components according to 
their weight scale, the highest of which are UUD config- 
urations. 
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